arXiv:1111.7030v2 [astro-ph.HE] 1 Dec 2011 


Astronomy & Astrophysics manuscript no. OFUpaperAA 
December 5, 201 1 


©ESO 2011 


Searching for soft relativistic jets in Core-collapse Supernovae 
with the IceCube Optical Follow-up Program 

R. Abbas®. y Abdoil®, T. Abu-Zayyad®, M. Ackerman r®. J. Adam.®, J. A. Aguilai®. M. Ahler®, M. M. Aller®, 
D. AltmanrP, K. Andeei®® J. Auffenbeq®, X. BaJ j3 l 44 l , M. Baket®, S. W. Barwick®, R. Bay®. J. L. Bazo Albd® 
K. Beattie®. J. J. BeatiyUHSl, S. Beche®, J. K. Becked, K.-H. Becked, M. L. Benabderrahmane®, S. BenZv®, 

J. Berdermanr®, R Berg ham®. D. Bciiey® E. Bcrnardin®. D. Bertram®, D. Z. Bessor®. D. B indig® M. BissokP, 
E. Blaufus®, J. BlumenthaP, D. J. BoersmaP. C. Bohn® D. Bose®, S. Bose®, O. Botne®, A. M. Browi®, 

S. Bukinl®, K. S. Caballero-Morr®. M. Carsoi®, D. Chirkir® B. Christy®, F. Clevermanr® S. Coher® 

C. Colnarc® D. F. Cowei®®, A. H. Cruz S i I vt®. ]yi . V. D’ AgostincP. M. Danninge®. J. DaughhctccP. J. C. Davi®. 

C. De Clercc®, T. Degne®, L. Demiror®, F. Descamp®, P. Desiat®, G. de Vries-Uiterweerc®, T. DeYomis®, 

J. C. Dfaz-Veled®, M. Dierckxsen®, J. Dreye®, J. P. Dumn®. M. Dunkmar®, J. Eiscl®, R. W. Ellsworth®, 

0. Engdegarc®, S. EulcP, P. A. Evensor®, O. Fadirar®, A. R. FazeljP, A. Fedyniteb®, J. Feintzeig® T. Feusel®, 

K. Filimono\P, C. Finley® T. Fischer- Wasel®, B. D. Fox® A. Franckowiak®*, R. Frankc® T. K. Gaissei® 

J. Gallaghe®, F. Gerhard® F. Gladstone®, T. GluserAami®, A. GoldschmidP, J. A. Goodmar®, D. Gort®, 

D. Gran®, T. Griese® A. Grot®®. S. G rul lor® M. Gurtne®, C. H®, A. Haj Ismai®, A. Hallgrei®, F. Halzer® 

K. Ha®, K. Hansoi®® D. HeineiP, K. Helbing® R. Hellaue®, P. Herque® S. Hickforc®, G. C. Hil®, 

K. D. Hoffmar®, B. HoffmanrP, A. Homeie® K. HoshinP, W. Huelsnid 18 ! 45 ! J.-P. Iliill.P. P. O. Hnltl®. 

K. HultqvisP, S. HussairP, A. Ishihar®, E. JacobP, J. JacobseiP, G. S. Japaridze^, H. JohanssoiP, 

K.-H. KamperP, A. Kappe®, T. Karg® A. Karlc® P. Ken in®. J. Kiryluk!® F. KislaP, S. R. Kleirl® 

J.-H. KohnP, G. KohnerP, H. KolanoskP, F. Kopk®, S. KoppeP, D. J. KoskinerP, M. KowalskP, T. Kowarik® 
M. Krasbei's®. G. KrolP, N. KurahashP, T. KuwabarP, M. Fabar®, K. FaihenP, H. FandsmarP, M. J. FarsorP, 

R. FaueiP, J. Fiinemanr® J. MadserP, A. MarottP, R. MaruyamP, K. MaseP, H. S. Mat®, K. MeagheP, 

M. MerckP, P. Meszarod 37 ! 38 ! T. Meure®, S. Miarecki® E. MiddelP, N. Milkd® J. MilleP, T. MontarulPMl 

R. MorseP, S. M. MoviP, R. NahnhaucP. J. W. NanP, U. NaumanrP, D. R. NygrerP, S. OdrowskP, A. OlivaP, 
M. Olive®, A. O’MurchadhP, S. PankniiP, L. PauP, C. Perez de los HeroP, J. PetrovicP A Pieg Sc® D. PielothP, 
R. PorranP. J. PosselP, P B. PriecP, G. T. PrzybylskP K. Rawlin.P P. RedP, E. Rescon®®.W. Rhode® 

M. Ribordy® M. RichmarP, J. P. RodrigueP, F. RothmaieP, C. RotP, T. RuheP, D. RutlcdgeP B. Ruzybayev® 
D. RyckbosclP, H.-G. SandeP, M. SantandeP, S. SarkaP, K. SchattcP, T. SchmidP, A. Schdnw akP, 

A. SchukrafP, A. SchulteP, O. Schul®® M. SchunckP, D. SeckeP, B. Semburj®, S. H. SecP, Y. SestaycP, 

S. Seunarine®, A. Silvestr®, G. M. SpiczakP, C. S pie ring® M. Stamatikosl 19 ! 49 ! T. Stance® T. StczelbcrgciP 
R. G. StokstacP, A. StossP, E. A. StrahleP, R. Stron®, M. Stiie®, G. W. SullivarP, Q. SwillenP, H. Taavol®, 

1. TaboadrP. A. TamburrcP, S. Tcr-AntonyaiP S. Tilav® P. A. ToalcP S. ToscancP, D. TosP, N. van EijndhoveiP, 
J. VandenbrouckeP, A. Van Overloot®, J. van SanterP, M. Vehrin® M. VogeP, C. WalckP, T. WaldenmaieP, 

M. WallrallP M. Waite®, Ch. Weave®, C. WendP, S. Westerhofl®, N. WhitehonP, K. WiebeP, C. H. WicbusclP 

D. R. William® R. WischnewskP, H. Wissing® M. WolP, T. R. WoocP, K. WoschnaggP C. XrP, D. F. XiP, 

X. W. XrP, J. P. Yanc/P G. YodlP. S. YoshidP, P. Zar/hitskyP, M. ZolP(the IceCube Collaboration), and 

C. W. AkerlolP S. B. Pandcy® F. YuarP, and W. Zheng® (the ROTSE Collaboration) 

(Affiliations can be found after the references) 

Preprint online version: December 5, 201 1 


ABSTRACT 

Context. Transient neutrino sources such as Gamma-Ray Bursts (GRBs) and Supernovae (SNe) are hypothesized to emit bursts of high-energy 
neutrinos on a time-scale of < 100 s. While GRB neutrinos would be produced in high relativistic jets, core-collapse SNe might host soft-relativistic 
jets, which become stalled in the outer layers of the progenitor star leading to an efficient production of high-energy neutrinos. 

Aims. To increase the sensitivity to these neutrinos and identify their sources, a low-threshold optical follow-up program for neutrino multiplets 
detected with the IceCube observatory has been implemented. 

Methods. If a neutrino multiplet, i.e. two or more neutrinos from the same direction within 100 s, is found by IceCube a trigger is sent to the 
Robotic Optical Transient Search Experiment, ROTSE. The 4 ROTSE telescopes immediately start an observation program of the corresponding 
region of the sky in order to detect an optical counterpart to the neutrino events. 

Results. No statistically significant excess in the rate of neutrino multiplets has been observed and furthermore no coincidence with an optical 
counterpart was found. 

Conclusions. The search allows, for the first time, to set stringent limits on current models predicting a high-energy neutrino flux from soft 
relativistic hadronic jets in core-collapse SNe. We conclude that a sub-population of SNe with typical Lorentz boost factor and jet energy of lip 
and 3 x 10 51 erg, respectively, does not exceed 4.2% at 90% confidence. 

Key words, neutrinos - supernovae - gamma-ray bursts 


1. Introduction 


High-energy astrophysical neutrinos are produced in interactions 
of charged cosmic rays with ambient photon or baryonic fields 


(for reviews se e |Anchordoqui & Montaruli|( 2010l; Chiarusi & 
Spurio] ( j2010j ); [Becker] (2008 i i; |Li pari | ( |20Q6 > ). Acceleration of 


these cosmic rays to very high energies takes place in astro- 
physical shocks. Neutrinos escape the acceleration region and 
propagate through space without interaction, while nuclei are 
deflected in magnetic fields and no longer point back to their 
source (for energies below ~10 2 ° eV). Unlike gamma-rays, neu- 
trinos are solely produced in hadronic processes and would 
therefore reveal the sources of charged cosmic rays. Gamma- 
ray bursts (see Meszaros (2002) and Zhang & Meszaros (2004) 


for reviews) could provide the environment and the required en- 
ergy to explain the production of the highest energy cosmic- 
rays (Waxmanj 1995). Recent observations indicate a connec- 
tion between long GRBs (duration > 2 s) and core-collapse 
supernovae (CCSNe). In several cases a gamma-ray burst or 
X-ray flash has been observed in coincidence with an optical 
SN light curve implying a common physical origin: a massive 
stellar explosion (see Woosley & Bloom) |2006 1 for a review). 
Furthermore, GRBs and CCSNe were found to release a compa- 
rable amount of kinetic energy. According to the collapsar model 


(MacFadyen & Woosley|1999[ Paczynski|1998| Woosley| 1993) 1, 
long GRBs have their origin in the collapse of a massive, rapidly 
rotating star into a black hole surrounded by an accretion disk. 
Relativistic jets with Lorentz boost factors of 100 - 1000 form 
along the stellar axis. The GRB-SN connection gives rise to the 
idea that GRBs and SNe might have the jet signature in com- 
mon and a certain fraction of core-collapse SNe might host soft 
relativistic jets with Lorentz boost factors around 5. The Lorentz 
boost factor of the jet might be determined by features of the pro- 
genitor star, such as its rotation. Compared to jets in GRBs, SN 
jets are suggested to be equally energetic but more baryon-rich, 
hence they are only mildly relativistic. Such soft relativistic jets 
would become stalled in the outer layers of the progenitor star, 
leading to essentially full absorption of the electromagnetic ra- 
diation emitted by the jet and, at the same time, an efficient pro- 
duction of high-energy neutrinos ( [Razza que et al.|2005| Ando & 
|Beaco m 2005 ). This motivates a search for neutrino emission, as 
neutrinos would be able to escape from within the star. 

The IceCube neutrino detector, located at the geographic 
South Pole, is built to detect high-energy astrophysical neutri- 
nos fAchterberg et al.|2006). So far GRB neutrino searches have 


been performed offline on AMANDA (Achterberg et al. 2008 
|2007bD and IceCube dAbbasi et all|2011b| |2009a| |2010) data. 


triggered by gamma-ray satellite detections. Time and direc- 
tion information provided by gamma-ray satellites allow an al- 
most background free search. An untriggered search was applied 
to AMANDA data (Achterberg et al. 2007b), which scanned 
the data for a clustering of neutrinos in time. Furthermore, a 
dedicated search for a neutrino signal in coincidence with the 
observed X-ray flash of SN 2008D has been conducted by 
IceCube (Abbasi et ak| 201 la) in order to test the soft jet sce- 
nario for CCSNe. Neither the GRB nor the SN neutrino search 
led to a detection yet, but set upper limits on the neutrino flux. 
Early SN detections, as in the case of SN 2008D, are very 
rare since X-ray telescopes have a limited field of view (FoV). 
However, neutrino telescopes cover half of the sky at any time. 
If neutrinos produced in soft relativistic SN jets are detected 
in real time, they can be used to trigger follow-up observa- 
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tions ( |Kowalski & Mohr||2007[ |Ageron et al.||20lT) . This is 
realized with the optical follow-up program presented here. 
Complementary to the offline searches, the optical follow-up 
program is an online search independent of satellite detections. 
It is sensitive to transient objects, including those which are ei- 
ther gamma-ray dark or not detected by gamma-ray satellites. In 
addition to a gain in sensitivity, the optical observations may al- 
low the identification of the transient neutrino source, be it a SN, 
GRB or any other transient phenomenon producing an optical 
signal. Hence, it enables us to test the plausible hypothesis of a 
soft relativistic SN jet and sheds light on the connection between 
GRBs, SNe and relativistic jets. 

In order to implement the optical follow-up program an online 
neutrino event selection was developed for IceCube. The data 
are processed online by a computer farm at the South Pole. A 
dedicated trigger selects neutrino burst candidates and the direc- 
tional information is transferred to the four ROTSE telescopes, 
which start the follow-up immediately and continue observations 
for several nights. The obtained optical data are searched for a 
transient counterpart. 

In the following, we present the optical follow-up program start- 
ing with section[2] which briefly describes the expected neutrino 
emission according to the soft SN jet model. Section [3] outlines 
the IceCube component of the program while section [4] focuses 
on the search for the optical counterpart. Finally, we discuss sys- 
tematic uncertainties in section [5] and show our results from the 
first year of data taking in section [6] with a focus on the SN soft 
jet model. We present a first limit on the hadronic jet production 
in CCSNe and conclude with a summary and outlook to future 
extensions of the program. 


2. SN neutrino flux 


Motivated by the GRB-SN connection Razzaque et al. ( |2005) 
proposed a model for high-energy neutrino production in soft 
relativistic CCSN jets. If protons are accelerated in the jets 
through Fermi acceleration in internal shocks, proton-proton- 
collisions will produce kaons and pions. The initial formula- 
tion of the model only considered neutrino production through 
pion decay (Razzaque et al. |2005) . It was extended by | Ando &| 
Beacom ( |2005) (hereafter AB05), who included neutrino pro- 
duction from kaon decay yielding a harder and hence more eas- 
ily detectable neutrino spectrum. AB05 present the calculation 
of the neutrino spectrum for a fixed Lorentz factor of To = 3 and 
a fixed jet energy of £j e t,o = 3 x 10 51 erg. In order to test a broader 
parameter space we calculate the neutrino flux as a function of 
the Lorentz boost factor I , the jet energy £jet and the density, p, 
of CCSN producing a jet. In the following, we derive the neu- 
trino flux for one SN at distance lOMpc assuming it hosts a jet 
pointing at us following the calculation of AB05. 

In the AB05 model, pions and kaons are produced with 20% 
of the parent proton energy and follow the original E 2 en- 
ergy spectrum originating from Fermi acceleration. However, 
protons lose energy through synchrotron radiation and inverse 
Compton scattering (radiative cooling) and through n p and Kp 
processes (hadronic cooling) causing a steepening of the spec- 
trum at higher energies. Above a certain break energy hadronic 
cooling becomes dominant and steepens the spectrum by a fac- 
tor E 1 while radiative cooling dominates above a second break 
energy resulting in a total suppression factor of E 2 . 

The daughter neutrino carries 25% of the pion energy or 50% of 
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the kaon energy and its energy is related to the meson energy in 
the jet framaMby 

E yM k) = TE' MK) / 4(2). (1) 

AB05 assume a hadronic cooling break of 30 GeV for pions and 
200 GeV for kaons. It depends strongly on the jet Lorentz factor 
T and the jet energy £j et : 

The radiative cooling break (0.1 GeV for pions and 20 GeV for 
kaons assumed by AB05) depends only on T: 

= ^0.1 (20) TeV. (3) 

l 0 

Note that for certain parameter configurations the order of 
the two break energies can change. Finally, the proton energy 
reaches its maximum at the photo-pion production threshold of 
E'p max = 7 x 10 4 GeV, where protons interact with the dense field 
of 4keV thermalized synchrotron photons. Neutrinos produced 
in resulting A + -decays are not considered by AB05. The cut-olf 
in the proton spectrum results in a cut-off in the neutrino spec- 
trum at 

^ 1 0.5(2 IjTeV. (4) 

1 0 

The normalization of the neutrino flux at 1 GeV scales with the 
jet energy. 

F vMK ), o = it 5 - (ft) 5 x 10 5 (5 x 10 1 2 ) GeV'W 2 , (5) 

£ jet.O U 0/ 

where the T 2 dependence is due to the assumed beaming with a 
jet opening angle of 6 oc 1 /T. Note that the expected number of 
neutrinos does not simply scale with the jet energy because the 
first break energy shifts with the jet energy. Depending on the 
choice of model parameters both the pion and kaon component 
of the spectrum can be hard (oc TV 2 ) or soft (oc TV 3 ) in the en- 
ergy range, which IceCube is sensitive to (TeV to PeV). Figure[I] 
illustrates the behavior of the total neutrino spectrum (sum of 
pion and kaon component) for different jet energies and Lorentz 
boost factors. 


3. IceCube 


The IceCube neutrino telescope has been under construction at 
the geographic South Pole since 2004 and was completed in the 
Antarctic summer of 2010/11. It is capable of detecting high en- 
ergy neutrinos with energies of <9(100) GeV and is most sen- 
sitive to muon neutrinos with energies in the TeV range and 
above. High-energy muon neutrinos undergoing charged current 
interactions in the ice or the underlying rock produce muons 
in neutrino-nucleon interactions. The muon travels in a direc- 
tion close to that of the neutrino and emits Cherenkov light, 
which is detected by a three dimensional array of light sensors. 
For this purpose, a volume of 1 km 3 of clear ice in depths be- 
tween 1450 and 2450 m below the ice surface was instrumented 
with 5160 digital optical modules (DOMs) attached to 86 ver- 
tical strings (|Achterberg et al. 2006). Each DOM consists of 


1 Variables in the comoving frame, i.e. the frame of the jet, are de- 

noted with a prime. 



Fig. 1. SN neutrino spectrum according to AB05 for one SN at 
distance lOMpc with the jet pointing to us. Different shades of 
gray indicate different Lorentz boost factors. Solid lines: £j et = 
3 x 10 51 erg. Dotted line: £j et = 0.3 x 10 51 erg. Dashed line: 
Tjet = 30 x 10 51 erg. 


a 25 cm diameter Hamamatsu photomultiplier tube (PMT) and 


supporting hardware inside a pressure glass sphere (Abbasi et al. 
|2009b| . The detector consists of 78 strings arranged in a hexag- 
onal shape with a string spacing of 125 m and DOMs separated 
vertically by 17 m, and 8 strings composing the low-energy ex- 
tension DeepCore ( Abbasi et al.|2011d] l, a densely spaced array 
in the bottom half of the detector. The observatory also includes 
a surface array, IceTop, to measure properties of extensive air 
showers and study the composition and spectrum of cosmic rays 
(Stane vp009) l. Figure [2] shows a top view of the IceCube detec- 
tor including DeepCore. Different colors/symbols indicate dif- 
ferent deployment stages. Before completion of the full detector, 
IceCube took data with the available number of strings. The opti- 
cal follow-up program has been fully operational since Dec. 16, 
2008. Here, we present the analysis of the data taken from Dec. 
16, 2008 to Dec. 31, 2009. Initially, 40 IceCube strings were tak- 
ing data (yellow upward-pointing triangle, green diamonds, red 
squares and magenta stars in Fig. [2|. In May 2009, additional 
19 strings were included (strings marked by purple downward- 
pointing triangles in Fig. [2] including the first DeepCore string 
marked by an open downward-pointing triangle). In the follow- 
ing, these deployment stages will be referred to as IC40 and 
IC59, respectively. 

To suppress the background caused by PMT noise or radioactive 
decay in the glass, a total number of eight DOMs with coinci- 
dent hits in a time window of 5 /-is are required for a trigger to be 
formed. A coincident hit is registered when a single DOM and 
its neighbor or next-to-neighbor DOM on the same string ex- 
ceed their charge threshold of 0.25 photoelectrons within a time 
window of 1 n s. If the detector is triggered the information of 
all triggered DOMs within a readout window starting 10 /rs be- 
fore the first hit and ending 10 /-is after the last hit is read out and 
merged to an event. Overlapping readout windows are merged 
together. The waveform of the PMT is digitized and sent to the 
surface. The waveforms have a length of up to 6.4 jus and can 
contain multiple hits. The total number of photoelectrons and 
their arrival times are extracted with an iterative unfolding algo- 
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Fig. 2. The IceCube detector: The full detector consists of 86 
strings with 60 DOMs attached to each string. Different col- 
ors/symbols indicate different deployment stages. The solid 
black line encircles the IC59 configuration, while the dashed line 
indicates the smaller IC40 configuration. 


rithm. 

The arrival time of the Cherenkov photons can be measured with 
an accuracy of ~3.3 ns (Abbasi et al .|2009b[ ). Several hit DOMs 
allow a reconstruction of the direction of muon-neutrinos with a 
precision of ~1°. To reduce the contribution from noise, only hits 
within a 6 /rs time window are used for the reconstruction (time 
window cleaning). This time window is defined as the window 
that contains the most hits during the event. Muons travel along a 
straight path through the detector while electrons or taus produce 
showers (cascades). Only the muon signature allows an accurate 
reconstruction of the direction. Directional information is cru- 
cial in order to provide coordinates for optical telescopes and we 
therefore consider only muon neutrino events. 


3. 1 . Online System 


In order to rapidly trigger optical telescopes, the first online 
analysis of high-energy neutrino events in IceCube was imple- 
mented. Unlike for the offline analyses, which are performed on 
an entire dataset (usually ~1 year of data) with time consuming 
reconstructions on a large computer cluster, the data are pro- 
cessed online by a computer cluster at the South Pole. During 
IC40 (IC59) data taking, the cluster consisted of approximately 
60 (100) processing nodes. The processing includes event re- 
construction and basic event selection. Several filters and recon- 


struction algorithms run on the cluster (see section 3.3 for details 


on the filters and reconstructions). During IC40 and IC59 events 
were read from the data acquisition and written to files with each 
file containing about 1 GB of data. The data volume of 1 GB cor- 
responds to a data taking period of -5 min or -4 x 10 5 events 
for IC40. Due to larger event sizes the same volume relates to 
-2.5 min or -2.5 x 10 5 events in IC59. The files were then dis- 
tributed to the processing nodes. Each file has to be processed by 
a single node. Although the processing time of ~35 ms per event 
is fast, the processing time per file amounts to ~4 h in IC40 and 
to ~2.5 h in IC59. For technical reasons, an additional latency of 
4 h was added to ensure correct ordering of the data after pro- 


cessing, which results in a total latency of 6.5 - 8 10 The pro- 
cessed and ordered data arrive on a dedicated machine (the anal- 
ysis client), which is not part of the parallel processing. There, 
a sophisticated event selection is applied based on the recon- 
structed event parameters. A multiplet trigger selects candidates 
of neutrino clusters, which are coincident in time and space. No 
further reconstruction algorithms need to be applied allowing a 
very fast filtering of the events. If a multiplet is found, its di- 
rectional information is transferred to Madison, Wisconsin, via 
the Iridium satellite network within < 10 s. From there, the mes- 
sage is forwarded to the four ROTSE telescopes via the Internet 
through a TCP-socket connection for immediate follow-up ob- 
servations. The stability and performance of the online system 
is constantly monitored in order to allow for a fast discovery 
of problems. To achieve this, test alerts are produced at a much 
higher rate (-100 test alerts per day compared to 25 real alerts per 
year) by the same pipeline and are also sent to the North. Their 
rate and delay time distributions are monitored using an auto- 
matically generated web page. If the rate deviates significantly 
from the expected, a notification is issued. 


3.2. Neutrino Data Set 


We present the analysis of data taken from Dec. 16, 2008 to Dec. 
31, 2009. This corresponds to a lifetime of 121 days with IC40 
and 186.4 days with IC59. Dead time is predominantly caused 
by calibration and commissioning runs during and after the con- 
struction season. The downtime of the online system amounts to 
6.8% mainly caused by downtime of the satellite transmission 
system. 

The background in a search for muon-neutrinos of astrophysical 
origin can be divided into two classes. One consists of atmo- 
spheric muons, created in meson decays in cosmic ray air show- 
ers, entering the detector from above. The other background is 
atmospheric neutrinos which originate in the same meson de- 
cays in cosmic ray air showers. Both are measured with IceCube 
and are well understood: the measurement of the atmospheric 
neutrino spectrum with IceCube in its 40-string configuration is 
discussed in |Abbasi et al. ( 201 lc| >, the atmospheric muon energy 
spectrum measured with the 22-string configuration is presented 
in Berghaus ( 2009] ). The flux of atmospheric muons exceeds the 
flux of atmospheric neutrinos by 5 orders of magnitude. The 
background of atmospheric muons can be reduced significantly 
by restricting the neutrino search to the northern hemisphere. 
Muons from the northern hemisphere cannot penetrate the Earth 
and reach the detector. However, a small fraction of the south- 
ern hemisphere muons are mis-reconstructed, i.e. truly down- 
going (entering the detector from above) but reconstructed as 
up-going. Owing to the large flux of atmospheric muons, these 
mis-reconstructed muons represent a significant contamination. 
Imposing requirements on the event reconstruction quality al- 
lows a suppression of the mis-reconstructed muon background. 


3.3. Neutrino Selection Criteria and Efficiency 

The expected neutrino signal according to the soft jet SN model 
can be calculated as a function of two model parameters: the 
Lorentz boost factor F and the jet energy £j et (see section |2V 
Signal events are simulated following the predicted neutrino mix 
spectrum in order to develop and optimize selection criteria to 
distinguish signal and background events. To suppress the back- 

2 With the start of operations with 79 strings the parallel processing 
was upgraded, reducing the total latency to a few minutes. 


4 


Abbasi et al.: IceCube Optical Follow-up Program 


ground of atmospheric neutrinos, which we cannot distinguish 
from the soft SN neutrino spectrum, we require the detection 
of at least two events within A t = 100 s and an angular differ- 
ence between their two reconstructed directions of AT < 4°. 
The choice of the size of the time window is motivated by the 
duration of the jet, i.e. the activity of the central engine, which is 
typically 10 s ( jRazzaque et al.|2005[ ). The observed gamma-ray 
emission from long GRBs has a typical length of 50 s ( |Gehrels 
et al. 2009| , which roughly corresponds to the time for a highly 
relativistic jet to penetrate the stellar envelope. The angular win- 
dow AT 1 is determined by the angular resolution of IceCube and 
was optimized along with the other selection parameters. The fi- 
nal set of selection cuts has been optimized in order to reach a 
background multiplet rate of ~25 per year corresponding to the 
maximal number of alerts accepted by ROTSE. Combining the 
neutrino measurement with the optical measurement allows the 
cuts to be relaxed compared to the neutrino point source analysis 
with IceCube (Abbasi et al .|201 le| yielding a larger background 
contamination but at the same time a higher signal passing rate. 
One doublet is not significant by itself, but may become signif- 
icant when the optical information is added. From the maximal 
allowed background multiplet rate a corresponding maximal sin- 
glet rate R \ can be estimated as follows. The probability to obtain 
a background triplet (three atmospheric neutrinos arriving by ac- 
cident within 100 s and within AT 1 ) or any multiplet of higher or- 
der is negligible, we therefore only consider doublets. Requiring 
no more than 25 background doublets per year (7\S < 25 year -1 ) 
corresponds to a rate of isotropic background singlets of: 


R\ < \ Ri 


Gnortb 

: At ■ AO 


7.2 mHz 
AT 1 ’ 


(6) 


paraboloid to the - log X space around the direction. Based on 
these more sophisticated reconstruction algorithms we select our 
final event stream. Level 3: a powerful parameter to reject mis- 
reconstructed events is given by the reduced negative logarithm 
of the likelihood, - log X/(Nch - 5) or a modified version given 
by - log -£/(Nch - 2) or - log £/(Nch - 3.5), where Nch is the 
number of triggered DOMs after time window cleaning. 

A large number of hits with a small time residual, N| )n -, i.e. 
registered within a time window [-25 ns, 75 ns] relative to the 
expected arrival time for unscattered light given by the track 
geometry, ensures a good track reconstruction quality, since 
photons causing those direct hits are less affected by scatter- 
ing. Furthermore, the maximum distance. Loir, along the recon- 
structed track direction between any two hits with small time 
residual ([-25 ns, 75 ns]) is a measure of the track quality. In ad- 
dition, only up-going events with zenith angle 6 > 85° (IC40) 
and 6 > 90° (IC59) are selectecQ In IC40, an additional cut on 
the doublet direction was applied selecting only doublets with 
a combined zenith of ^Doublet > 90°. A single-iteration likeli- 
hood reconstruction (llhl) and a second two-iteration likelihood 
reconstruction (llh2) using the seed track and the inverted seed 
track of llhl were used in IC40. The ten-fold iterative likelihood 
reconstruction (lOit) was used in both IC40 and IC59. 

The final cuts at Level 3 for IC40 are 

6>n h i > 85° AND (? llh2 > 85° AND 0 1Oll > 85° AND 
- log£/(N C h - 5) < 8.85 AND 
((N Dir > 7 AND L Dir > 225) OR N ch > 200) 

and for IC59 

6> I0i , > 90° AND - log -£/(Nch - 3.5) < 7.7 AND 
((N Dir > 7 AND L Dir > 250) OR N ch > 100). 


where AQ = tt(A'V) 2 is the solid angle defined by the dou- 
blet condition and Qr>jorth = 20627(°) 2 is the solid angle of the 
Northern sky, i.e. 27r. The event selection is optimized in or- 
der to restrict the singlet rate to 7.2 mHz/ AT while obtaining 
the best signal passing rate. AT = 4° was found to be the best 
choice during the cut optimization. This corresponds to a sin- 
glet rate of 1.8 mHz. The event selection is based on one of 
the standard IceCube muon event filters (in the following re- 
ferred to as Level 1), which is commonly used by several of- 
fline analysis. It is discussed in detail in Abbasi et al. ( 2011e| . 
The Level 1 filter selects muon tracks, which are reconstructed 
as up-going (passing through the Earth) based on fast and simple 
algorithms. It selects ~2% of all triggers with a signal efficiency 
of 90% for up-going neutrinos following an E 2 -spectrum. It is 
still largely dominated by atmospheric muons. To the Level 1 
filter stream we apply additional selection criteria mainly based 
on the reconstruction quality. This yields the so-called Level 2 
filter stream. The significantly smaller rate of -3 Hz allows us to 
perform time consuming reconstructions online, which provide 
a more accurate estimate of the event direction. These recon- 
structions are based on a muon-likelihood function described in 
Ahrens et al. j ( 2004) >, which parametrizes the probability of ob- 


serving the spatial distribution and timing of the hits with re- 
spect to a muon track hypothesis. The negative logarithm of 
the likelihood, -logX, is minimized, i.e. the likelihood, X, is 
maximized by varying the track direction to yield the best-fit 
direction and position for the muon track. Iterative fits repeat 
the minimization with a different initial track hypothesis to re- 
duce the problem of local minima. The iteration with the small- 
est minimum is the final fit result. At Level 2, a ten-fold it- 
erative likelihood reconstruction is available. The uncertainty 
on the reconstructed direction, cr, is obtained from a fit of a 


The spectrum obtained from the AB05 model can be either hard 
or soft depending on the choice of model parameters. The choice 
of cuts specified above is a compromise yielding adequate pass- 
ing rates for all considered model parameters. Figure [3] shows 
the energy dependence of the Level 3 filter efficiency relative to 
Level 2 for well-reconstructed events (|~'? lrul , - 4^ rcc0 [ < 3° and 
$true > 90°, where the unit vector^? indicates the track direction). 
The filter efficiency is defined as the fraction of simulated signal 
events passing the filter. For energies above 100 TeV the filter is 
90% efficient while for lower energies the efficiency decreases to 
50% at 1 TeV and 20% at 100 GeV. The cuts were adjusted from 
IC40 to IC59 to account for the larger detector volume in order 
to obtain a similar data passing rate. Table[l]shows the data pass- 
ing rates for IC40 and IC59 at different cut levels. Furthermore, 
the table contains the expected number of detected SN neutri- 
nos for a SN at distance c/sn = lOMpc with a jet of energy 
£jet = 3 x 10 51 erg, which points towards Earth, and two choices 
of the boost factor L (4 and 10). The expected number of well- 
reconstructed SN neutrinos is given in brackets since only these 
events are useful to trigger optical telescopes. Figure[4]shows the 
expected number of events for different model parameter config- 
urations for IC59. Table[l]shows that we expect many more neu- 
trino events for a SN with a high boost factor. On the other hand, 
a large boost factor implies a small jet opening angle (6 oc 1 /F) 
and hence a smaller probability of the jet pointing towards Earth. 
Furthermore, one sees an increase in the expected number of SN 
neutrinos at Level 3 from IC40 to IC59. The detector volume 
increased by roughly 50%. However, owing to improved perfor- 
mance of the reconstruction algorithms applied to data of the 

3 Southern hemisphere events have zenith angles of 0° < 6 < 90° in 
the coordinate system of IceCube. 
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Table 1. Cut Levels and Event Rates. Values in brackets refer to well-reconstructed events. 


Cut Level 

Event Rate 

Neutrino events for SN at d s n = 

lOMpc 




r = 4, £ jet = 3 x 10 51 erg 

T = 10, £j e t 

= 3 x 10 51 erg 


IC40 

IC59 

IC40 IC59 

IC40 

IC59 

Level 1 

20.7 Hz 

22.7 Hz 

68.0(18.6) 133.3 (34.7) 

3385.5 (1081.3) 

5304.0 (1877.7) 

Level 2 

2.74 Hz 

3.32Hz 

48.1 (17.6) 100.0(32.9) 

2544.6 (1076.0) 

4225.9 (1801.1) 

Level 3 

2. 17 mHz 

1.86 mHz 

13.3(8.7) 22.4(16.3) 

947.4 (674.6) 

1441.7 (1153.0) 



Fig. 3. Filter efficiency of Level 3 relative to Level 2 as function 
of energy for well-reconstructed events (IC59). 
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Fig. 5. Neutrino effective area at the final selection level for well- 
reconstructed events for IC40 (blue squares) and IC59 (red cir- 
cles). 
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Fig. 4. Expected number of SN neutrinos from a SN at distance 
lOMpc with a jet pointing towards Earth for different model pa- 
rameter combinations (IC59). 


larger 1C59 detector, the signal passing rate at Level 3 increased 
by 52% (E = 10) and 68% (r = 4). The increase for well- 
reconstructed events is even higher (71% for T = 10 and 87% 
for T = 4). The neutrino effective area for well-reconstructed 
events at final cut level is shown in Fig. [5] for IC40 and IC59. 
The Level 3 cuts reduce the singlet rate to ~1.8 mHz for IC59. 
The IC40 rate is slightly higher at 2.2 mHz, but an additional cut 
on the doublet direction (^doublet < 90°) ensures a multiplet rate 
of less than 25 per year. The Level 3 data stream consists of 37% 
(70%) atmospheric neutrinos for IC40 (IC59), the rest is a con- 
tamination of mis-reconstructed atmospheric muons. In addition 



Fig. 6. Angular difference between the reconstructed directions 
of two neutrinos with identical true direction for jet parameters 
r = 3 and E jct = 3x 10 51 erg (IC59). 


- as described above - we require at least two events to arrive 
within At = 100 s and with an angular distance of AT < 4° to 
reduce the background of atmospheric neutrinos. The signal ef- 
ficiency of the angular coincidence cut for different zenith bands 
is displayed in Fig. [6] The average passing rates depend on the 
assumed energy spectrum (i.e. the model parameters) and are in 
the range of 56-69% (IC40) and 60-74% (IC59) for 1 < T < 10 
and 3.1 x 10 49 erg < E; et < 3 x 10 53 erg. The reconstruction ac- 
curacy has improved because of the increased detector volume. 
Large reconstruction uncertainties might lead to mis-pointing 
of the telescopes and in the worst case the real source position 
might lie outside ROTSE’s field of view of 1.85° x 1.85°. To im- 
prove the accuracy of the direction forwarded to the telescopes. 
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Fig. 7. Doublet resolution using an ordinary mean (blue circles), 
a weighted mean (red squares) compared to the singlet resolution 
(green triangles) for signal neutrinos (r = 3, £jet = 3 x 10 51 erg). 
Applying the directional coincidence cut A'F < 4° (black tri- 
angles) keeps mainly well-reconstructed doublets and yields a 
further improvement of the doublet resolution (IC59). 



Fig. 9. Efficiency to find inserted SNe as a function of the appar- 
ent peak magnitude averaged over all alert directions assuming 
a template light curve. 
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Fig. 8. Deviation of reconstructed multiplet direction from tme 
direction (IC59). The black box shows ROTSE’s field of view of 
1.85° x 1.85°. 


the doublet direction is calculated as a weighted mean from the 
directions of the individual events in the multiplet. The single 
events are weighted with 1/cr 2 , where cr is the reconstruction er- 
ror estimated by the paraboloid fit. Compared to single events, 
doublets have a better resolution. The weighting improves the 
doublet resolution as can be seen in Fig. [7] It further improves 
after applying the directional coincidence condition AT* < 4°. 
Figure [8] shows the doublet point spread function in 0-<p - space 
compared to ROTSE’s FoV. Depending on the model parame- 
ters, 41-53% (IC40) and 44-61% (IC59) of all doublet events lie 
within ROTSE’s field of view, providing a good match for this 
search. 


4. Search for Optical Counterparts 

The IceCube multiplet alerts are forwarded to the robotic opti- 
cal transient search experiment (ROTSE), which consists of four 
identical telescopes located in Australia, Texas, Namibia and 
Turkey (Akerlof et al. 2003} >. The telescopes stand out because 
of their arge FoV of 1.85° x 1.85° and a rapid response with a 
typical telescope slew time of 4 s to move the telescope from the 


standby to the desired position. The telescopes have a parabolic 
primary mirror with a diameter of 45 cm. To be sensitive to weak 
sources no bandwidth filter is used. ROTSE is most sensitive in 
the R-band (-650 nm). The wide field of view is imaged onto a 
back-illuminated thinned CCD with 2048 x2048 13.5 /vrn pixels. 
The camera has a fast readout cycle of 6 s and is cooled to -40° C 
in order to reduce thermal noise. For a 60 s exposure at optimal 
conditions the limiting magnitude is around Mr ~ 18.5, which is 
well suited for a study of GRB afterglows during the first hour 
(up to one day for very bright afterglows) and SN light curves 
with apparent peak magnitude <16. The corresponding FWHM 
(full width at half maximum) of the stellar images is < 2.5 pix- 
els (8.1 arcseconds). The telescopes are operated robotically and 
managed by a fully-automated system. Observations are sched- 
uled in a queue and are processed in the order of their assigned 
priority. IceCube triggers have second highest priority after GRB 
follow-ups triggered by the GRB Coordinate Network (GCN). 
Once an IceCube alert is received by one of the telescopes, a 
predefined observation program is started and the correspond- 
ing region of the night sky will be observed within seconds: 
the prompt observation includes thirty exposures of 60 seconds 
length, follow-up observations are performed for 14 nights. This 
was extended on Oct. 27, 2009 to 24 nights, with daily obser- 
vations for 12 nights and then observations during every second 
night up to day 24 after the trigger was received. Eight images 
with 60 seconds exposure time are taken per night. The prompt 
observation is motivated by the typical rapidly decaying light 
curve of a GRB afterglow, while the follow-up observations dur- 
ing 14 or 24 nights permit the identification of a rising SN light 
curve. With IC40 and IC59, the online processing latency of sev- 
eral hours made the search for an optical GRB afterglow unfea- 
sible. We therefore focus on the SN light curve detection in the 
ROTSE data. 

Image correction and calibration are performed at the telescope 
sites. The images of each night are co-added in order to obtain a 
deeper image. Co-adding includes a geometrical transformation 
to correct for optical distortions and mis-alignment of the images 
before the pixel contents are added. A reference image is sub- 
tracted from each co-added image. As deep images are usually 
not available for the positions we would like to observe, we ini- 
tially choose the deepest image of our observing sequence as the 
reference image (due to different weather conditions in different 
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nights the image quality varies). The source could be present in 
the reference image, but it would be of different brightness com- 
pared to earlier or later images, and both positive and negative 
deviating sources in the subtracted images will be detected. If 
no early image of good quality was available to be selected as 
reference image (30% of the alerts) we take another deep image 
several months later. Both SN light curves and GRB afterglows 
would have faded after a few weeks and would not be present in 
the newly taken reference image. 


The applied subtraction algorithm was developed by Yuan & 
Akerlof| ( 2008| >. Both the new and the reference image are folded 
by a kernel function in order to match their point spread func- 
tions (PSFs). The convoluted images then allow a pixel by 
pixel subtraction. The software SExtractor (source extractor), 
Bertin & Arnouts] ( 1996 1, extracts all objects from an image by 
first determining the background and then identifying clusters 
of pixels with a significance > let above the background level. 
All extracted objects found in the subtracted images with signal- 
to-noise ratio larger than 5 are candidates for variable sources. 
However, bad image quality, failed image convolution, bad pix- 
els and other effects frequently cause artifacts in the subtraction 
process, requiring further selection of the candidates. A candi- 
date identification algorithm including a boosted decision tree 
is applied to classify candidates according to geometrical and 
variability criteria. The algorithm was trained using a signal of 
a SN light curve starting at the time of neutrino detection. We 
use a SNIbc template light curve based on SN1999ex (Hamuy 
et al.[2002$ . The time of the shock break-out was measured for 
SN1999ex and provides a time stamp for the explosion, i.e. the 
start of the light curve, which can be associated with high-energy 
neutrino emission. To simulate the SN, fake stars are inserted on 
top of galaxies in every single image from the observation se- 
quence with a brightness according to the SN light curve tem- 
plate. To ensure that the PSF of the inserted star reflects all the 
features of the PSF of existing stars in the image we use the 
mean PSF of all stars in a 291 x 291 pixel box around the in- 
sertion coordinates. The PSF of a single star is obtained from a 
box of 15 x 15 pixel around the star’s center. The manipulated 
images are processed with the same pipeline as described above. 
The signal sample consisting of inserted fake SNe is used to train 
the classification algorithm as well as to estimate ROTSE’s ef- 
ficiency. The efficiency is given by the fraction of inserted SNe, 
that has been detected by the processing and candidate identi- 
fication. For some inserted SNe the detection algorithm fails: 
If the quality of the image is bad (e.g. large average FWHM 
or small limiting magnitude) the image convolution performs 
badly. Candidates close to saturated objects or close to objects 
listed in the two micron all sky survey (2MASS) point source 
catalog, which are very likely stars, are removed automatically. 
Figure [9] shows the fraction of simulated SNe that are found by 
the algorithm as a function of the inserted apparent peak magni- 
tude. The efficiency as a function of the apparent peak magnitude 
can be converted to the efficiency as a function of SN distance 
srotse(^) (see Fig.[l0|) assuming an absolute R-band magnitude 
of M = -18 ± 1 for core-collapse SNe ( jRichardson et al. 12006) 1. 
The relation of distance and magnitude is given by 


(7) 


m = M + 5 ( log 10 — - 1 


where m is the relative magnitude. 

The efficiency is used to calculate the number of expected SNe 
detections for a given signal neutrino hypothesis. However, it 
may also be used to estimate, for a given SN rate, the number 
of accidental coincidences between a neutrino doublet and a SN 



Fig. 10. Black curve: Efficiency to detect core-collapse SN as 
a function of the distance to the SN assuming an absolute R- 
band magnitude of -18 ± 1. Shaded region: lower bound assum- 
ing an absolute magnitude of -17, upper bound assuming -19. 
The breaks in the shaded regions are connected to the binning in 
Fig . [9] The binning effect is washed out in the black curve due to 
the assumed uncertainty in the absolute magnitude distribution 
of+1. 


detectable by ROTSE (see section [6]). After application of the 
classification algorithm -10 - 200 candidates for variable ob- 
jects remain per alert depending on the quality of the images 
and the Galactic latitude. Fields close to the Galactic plane con- 
tain a large number of stars, which complicates the subtraction 
and thus results in more candidates caused by subtraction arti- 
facts. Tightening the cuts would reduce the number of candidates 
but at the same time reduce our sensitivity. The final candidates 
are summarized on a web page and are inspected visually by 
people trained on the signal simulation. The web page shows 
a 100 x 100 (pixel) 2 extract of the new, the reference and the 
subtracted image for each night. For comparison, an image from 
the digitized sky survey (DS^]) of the same patch of the sky 
is shown, which is deeper than the ROTSE images. In addition 
links to cat alogs, such as NEEp] 2MASS ( Skrutskie et al. 12006 


and SDSS (Adelman-McCarthy et al.|2008p are provided. On the 
basis of this information the scanners have to decide whether the 
candidate is a SN, a variable star or a subtraction artifact. SN 
candidate identification by the human eye works well as shown 
in the galaxy zoo SN project ( (Smith et ak 2010 1 . The visual scan- 
ning of the final candidates was carried out by three individual 
persons, who obtained similar results. All simulated SNe, which 
passed the computer selection, were identified in the scanning, 
i.e. the efficiency was 100%. Also the rate of false positives is ex- 
pected to be small, because for a potential SN candidate its light 
curve and host galaxy properties would be inspected in detail. 
We hence assume, that the systematic uncertainty introduced by 
the scanning process is negligible. Note that in the future, for 
candidates identified in real time a spectrum can be obtained to 
ensure an unambiguous identification of the SN. 


5. Systematic Uncertainties 

Both the simulated neutrino sensitivity and the SN sensitivity are 
subject to systematic uncertainties. These systematic uncertain- 
ties are included in the limit calculation. In this limit calculation, 

4 

5 


http : // archive . stsci . edu/dss/ 
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Fig. 11. Systematic uncertainties (relative to the predicted SN 
neutrino event rate) depending on model parameters F and £j et . 


Monte Carlo experiments are performed drawing the number of 
signal neutrino-events following a Poisson distribution (see ap- 
pendix [A] for details on the limit calculation). Systematic un- 
certainties are included by smearing the Poisson mean, i.e. the 
Poisson mean is multiplied by a factor following a Gaussian dis- 
tribution with mean one and a width given by the systematic 
uncertainty. 


Table 2. Systematic Uncertainties and their Influence on the 
Event Rates. 


Source of 
uncertainty 

Uncertainty 

Average influence 
on event rate 

DOM Efficiency 

±10% 

±18% 

Photon propagation 

±10% 

±13% 

Neutrino cross section 

[-3%, -10%] 

-6% 

Muon energy loss 

±1% 

±1% 


Table 3. Systematic Uncertainties 


Source of 
uncertainty 

Uncertainty 

Influence 
on SN rate 

Magnitude 

measurement 

0.3mag 

-17.6%, +26.6% 

CC SN rate 

±30% 

±30% 

Inhomogeneity of 
the local universe 

±30% 

±30% 

Quadratic sum 


-45.9, +50.1% 


5.1. Systematic Uncertainties on the Neutrino Event Rate 


The main sources of systematic uncertainties arise in the de- 
scription of the DOM efficiency and the photon propagation in 
ice. The systematic uncertainties due to photon propagation are 
evaluated by performing dedicated simulations with scattering 
and absorption coefficients varied within their uncertainties of 
-10% (Ackermann et al.][2006| l. The maximum difference was 
found between the case where both scattering and absorption 
were increased by 10% and the case where both were decreased 
by 10%. Varying the DOM efficiency resulted in a variation of 
the event rate of up to 18%. The neutrino cross section used 
in the neutrino simulation is based on CTEQ5 measurements, 
which are not up-to-date anymore. The latest cross section cal- 
culation by |Cooper-Sarkar & Sarkar ( 2008 1 differ from the cross 
sections used in the IceCube neutrino simulation by up to 10% 
in the relevant energy regime. To first order, the neutrino rate 
depends linearly on the cross section. Folding the expected SN 
neutrino spectrum with an energy dependent correction factor 
allows us to calculate the effect on the neutrino event rate, which 
is up to 6%. 

Finally, the uncertainty in the muon energy loss amounts to 1%, 
resulting in a 1% influence on the event rate (Achterberg et al. 
2007 a) l. The systematic uncertainties are energy dependent and 
the SN neutrino spectrum varies with the model parameters. 
Therefore, we calculate the quadratic sum of all listed system- 
atic uncertainties for each combination of model parameters (see 
Fig. 111. Table [2] summarizes the systematic uncertainties con- 
sidered in this analysis and their influence on the event rate for 
those model parameters, where the effect is largest. 


5.2. Systematic Uncertainties on the SN Sensitivity 

The number of expected observed SNe depends on the sen- 
sitivity of the telescopes. The estimate described in section [4] 
yields the efficiency as a function of the SN peak magnitude. 


The photometric zero-points as determined from USNO A2.0 11- 
band magnitudes have typical systematic uncertainties of up to 
0.30 mag (Rykoff et al. ( 2005| ) and references therein). Shifting 
the efficiency curve by +0.3 mag results in a variation of the 
expected number of SNe measured by ROTSE of [-17.6%, 
+26.6%]. 

The expected number of accidentally found SNe depends on 
the overall core-collapse SN rate, which is assumed to be 1 
SN per year within a sphere of radius lOMpc (continuum limit 
from Ando et al. ( 2005| >). The true SN rate might be higher since 
nearby SNe surveys tend not to target small galaxies. During the 


last decade 17 SNe within lOMpc were observed (Kistler et al. 


2008 1. We assume a systematic uncertainty of 30% due to in- 


homogeneity of the local universe and 30% on the CCSN rate. 
Systematic uncertainty introduced by the scanning process are 
considered negligible. 


6. Results 

This paper presents the results from the analysis of data taking in 
the period of Dec. 16, 2008 to Dec. 31, 2009. IceCube was run- 
ning initially in the IC40 configuration (Dec. 16, 2008 to May 
20, 2009) and later in the IC59 configuration (May 20, 2009 to 
Dec. 31, 2009^ Table [ 4 ] shows the number of detected and ex- 
pected doublets and triplets for the IC40 and the IC59 datasets 
as well as the number of detected and expected optical SN coun- 
terparts. The IceCube expectation based on a background-only 
hypothesis was obtained from scrambled datasets. To correctly 
incorporate detector asymmetries, seasonal variations and up- 
time gaps, we used the entire IC40 and IC59 datasets and ex- 


6 Note that the IceCube detector was running in the 40-string con- 
figuration already before Dec. 2008 and took data with the 59-string 
configuration also after Dec. 2009. The 2010 dataset is currently being 
analyzed. 
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Table 4. Measured and Expected Number of Multiplets 



SN 

Doublets 

Triplets 



IC40 

IC59 

IC40 

IC59 

measured 

0 

15 

19 

0 

0 

expected 

0.074 

8.55 

15.66 

0.0028 

0.0040 


changed the event directions randomly while keeping the event 
times fixed. For each scrambled dataset we obtain the number of 
doublets by comparing event directions. The number of doublets 
in both data sets (IC40 and IC59) shows a small excess, which 
corresponds to a 2.1 cr effect and is thus not statistically signif- 
icant. To estimate the expected number of randomly coincident 
SN detections, we assume a core-collapse SN rate of 1 per year 


within a sphere of radius lOMpc, i.e. 2.4x 1CF 4 y _1 Mpc , and a 
Gaussian absolute magnitude distribution with mean of - 1 8 mag 
and standard deviation of 1 mag (Richardson et al.|2006| >. Based 
on the efficiency estimated in section [4] we can calculate the 
rate of core-collapse SNe that could be detected by ROTSE, if 
it would continuously survey the full sky, = 3823 y 

For this we integrated the CCSN rate over the accessible volume 


weighted with the efficiency displayed in Fig. 10 The number 


of expected accidental SN detections (i.e. a SN detection in co- 
incidence with a background neutrino multiplet) is 


-Nsn ,exp — • 2V alerts 


n 


ROTSE 


a 


r» 

• K 


■sky 


ROTSE 

CCSN 


= 0.074 


( 8 ) 


where /V a | erts = 17 is the number of multiplet alerts followed-up 
by ROTSE, Grotse = 1.85° x 1.85° is the solid angle covered 
by ROTSE’s field of view and G s k y = 41253(°) 2 is the all sky 
solid angle. ATsn is the time window in which we accept a co- 
incidence of neutrino and optical signals. It has to be larger than 
the uncertainty of the SN explosion time. In Cowen et al.]( |2010| > 
it is shown, that the explosion time can be estimated with an 
accuracy of ~ 1 day if early data are available. We choose to be 
conservative and take A7 'sn = 5 days. In total, 31 alerts were 
forwarded to the ROTSE telescopes. Five could not be observed 
because they were too close to the Sun. For two alerts no good 
data could be collected. Seven alerts were discarded because the 
corresponding fields were too close to the Galactic plane and 
hence too crowded. Thus 17 good optical datasets remained for 
the analysis. The data were processed as described above. No 
optical SN counterpart was found in the data. We calculate the 
limit on the AB05 model parameters following the description in 
appendix |A] for the jet Lorentz boost factors T = 6, 8, 10 and in 
each case vary the jet energy Ejet and the rate of SNe with jets p. 
The algorithm was formulated prior to the start of the program. 
The systematic uncertainties discussed in section [5] are included 
in the limit calculation. For each I -value the 90% confidence 
region in the £j ct -p-plane is displayed in Fig. [12] The colored re- 
gions are excluded with 90% confidence. The limits include the 
optical information, i.e. that no optical counterpart was found. 
This improved the limit and allows tests of 5-25% smaller CCSN 
rates. The largest improvement is obtained for small jet energies 
and large CCSN rates. The most stringent limit can be set for 
high Lorentz factors, while for small F the constraints are weak. 
At 90% confidence level, a sub-population of SNe with typical 
values of T and Ej e t of 10 and 3 x 10 51 erg, respectively, does not 
exceed 4.2%. This is the first limit on CCSN jets using neutrino 
information. 



IOg io( E jet[ erg ]) 


Fig. 12. Limits on the choked jet SN model Ando & Beacomi 
( j2005| for different Lorentz boost factors T as a function of the 
rate of SNe with jets p and the jet energy £j et . The colored re- 
gions are excluded at 90% confidence level. Horizontal dashed 
lines indicate a fraction of SNe with jets of 100%, 10% or 1% 
(relative to an assumed CCSN rate of 1 per year within a sphere 
of radius 10 Mpc). 


7. Summary and Outlook 

This first analysis, using the four ROTSE telescopes, proves the 
feasibility of the program for follow-up observations triggered 
by neutrino multiplets detected by IceCube. The technical chal- 
lenge of analyzing neutrino data in real time at the remote lo- 
cation of the South Pole and triggering optical telescopes has 
been solved. First meaningful limits to the SN slow-jet hypoth- 
esis could be derived already after the first year of operation. 
Especially in cases of high Lorentz boost factors of T = 10 strin- 
gent limits on the soft jet SN model are obtained. |Soderberg et al.l 
( 20 1 0j > obtain an estimate on the fraction of SNe harboring a 
central engine, which powers a relativistic outflow, from a ra- 
dio survey of type Ibc SNe. They conclude that the rate is about 
1%, consistent with the inferred rate of nearby GRBs. Our ap- 
proach is completely independent and for the first time directly 
tests hadronic acceleration in CCSNe, while the radio counter- 
part is sensitive to leptonic acceleration. 

The instrumented volume of IceCube has now increased to a 
cubic kilometer yielding an increased sensitivity to high-energy 
neutrinos. In addition the live time is growing continuously. The 
delay of processing neutrino data at the South Pole has been re- 
duced significantly from several hours to a few minutes. This 
results in the possibility of a very fast follow-up and allows the 
detection of GRB afterglows, which fade rapidly below the tele- 
scope’s detection threshold. 

In addition, a single high-energy neutrino event trigger (in ad- 
dition to the mutliplet trigger) is under development, which will 
further increase the sensitivity of the program especially for hard 
GRB neutrino spectra. 

Because of the successful operation of the optical follow-up pro- 
gram with ROTSE, the program was extended in August 2010 to 
the Palomar Transient Factory (PTF) (Law et al.|2009 ; Rau et ah 
2009 ), which will provide deeper images and a fast processing 
pipeline including a spectroscopic follow-up of interesting SN 
candidates. Furthermore, an X-ray follow-up by the Swift satel- 
lite ( jGehrels et al.||2004 1 of the most significant multiplets has 
been set up and started operations in February 2011. 
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Appendix A: Limit Calculation 

Motivated by the GRB-SNe connection, the soft SN jet model 
predicts the production of mildly relativistic baryon loaded jets 
in core-collapse SNe. The resulting neutrino flux depends on the 
jet energy £j et and the jet boost factor T. The rate of detectable 
SNe depends on the rate of core-collapse SNe producing a jet p. 
In order to test the model we define a test statistic T consisting 
of an IceCube term Tic and a ROTSE term Trotse- 
The probability to detect A'; ; or more multiplet events in IceCube 
with multiplicity k over a background expectation of p/ : is given 
by the sum over Poisson probabilities: 


P(N k ,p k ) = J] T e " W - 

i=N k l ' 


(A.l) 


Combining all multiplicities and the two datasets yields the test 
statistic Tic 

CO 

T IC = Y] P(Vf 40 ,/4 C40 ) • P(Nl c 59 ,pl c59 ). (A. 2) 

k = 2 

In addition to the IceCube information (i.e. number of doublets 
and multiplets of higher order) we add information obtained by 
the optical observations. The probability to observe N$n or more 
optical SN counterparts based on the expected number p SN of 
accidentally observed SN in coincidence with an IceCube multi- 
plet (given by Eq.[8]) is given by the sum of Poisson probabilities: 

P SN (N SN ,p SN ) = J] (A. 3) 

i = n sn *' 


If one or more optical counterparts were observed the signifi- 
cance could be improved by adding neutrino timing information 
as well as the distance information of the object found. The prob- 
ability P, to find a time difference of At or smaller between the 
first and the last neutrino event in the multiplet (defined by the 
maximal temporal difference of 100 s between the neutrino ar- 
rival times) due to a background fluctuation assuming a uniform 
background is given by 


P,= 


At 

100 s' 


(A. 4) 


Hence, assuming a generic SN prediction of a 10 s wide neutrino 
pulse results in a factor of 10 lower chance probability. 

Taking into account the SN distance allows us to compute the 
probability Pd to observe a background SN at a distance d < cIsn 


Pd 


P'c?lf(d < d SN ) 

dROTSE 

“CCSN 


(A. 5) 


where < dsn) is the rate of SN observable by the 

ROTSE telescopes within a sphere of radius dsn ■ ^ccsn E ' s t ^ le 
total number of SN observable by ROTSE. Accidental coinci- 
dences will be distributed following the SN rate (i.e. the square 
of the distance) folded with ROTSE’s sensitivity as a func- 
tion of distance srotseW) (see section |4j. Signal events have 
a strong preference to close-by SNe, since only these will lead 
to a neutrino flux large enough to produce a detectable multiplet 
in IceCube. While ROTSE can only detect close-by SNe, more 
powerful telescopes can access a much larger volume and would 
essentially always detect a SN in their field of view. Hence, 
the additional factor Pd becomes important to account for the 
SN distance. The additional terms P, and P,i for each observed 
SN light curve are combined with Psn yielding the test statistic 
Trotse- 


Nsn 

Trotse = Psn(N SN , p SN ) f| P u ■ P di . (A. 6) 

i=i 

Combining all information into one test statistic T yields: 

T = Tic • Trotse- (A. 7) 

To obtain a proper confidence region for exclusion of the 
model we perform 10000 Monte Carlo (MC) experiments for 
each combination of model parameters. The model prediction, 
V ( f ,/IC59 and (Vf N , depends on the model parameters, £ jct , T 
and p. is obtained from the neutrino signal simulation 

weighted with the corresponding AB05 spectrum. The neutrino 
spectrum varies with Ejet and T as presented in section [2] The 
number of predicted SNe, /V 9N , depends on the number of neu- 
trino multiplets, i.e. number of telescope pointings, folded with 
the sensitivity of the telescope. In our signal estimation we have 
not accounted for mixed multiplets due to a single SN neutrino 
in coincidence with a background neutrino. While in principle, 
these can be indentified through an optical counterpart, we esti- 
mate the rate to be at most a few percent of the pure, signal-only 
multiplets. We hence neglected the extra contribution. The aver- 
age number of background multiplets (V[ ( j 40 ^ IC59 is known from 
scrambling. 

For each MC experiment the number of signal and background 
multiplets are drawn following a Poisson distribution. In case of 
the signal, the systematic uncertainties are included by smearing 
the Poisson mean, i.e. the Poisson mean is multiplied by a fac- 
tor following a Gaussian distribution with mean one and a width 
given by the systematic uncertainties summarized in section [5] 

If an optical counterpart is drawn in the MC simulation (Nf"> 
1) we calculate the additional terms P t and P,i following equa- 
tions IA.4I and IA.5I The time difference between the SN neutri- 
nos is set to At = 10 s and the SN distance is thrown following a 
spatially isotropic distribution folded with ROTSE’s efficiency. 
For each MC experiment T is calculated following equation |A.7| 
The fraction of MC experiments resulting in a smaller value of T 
than that of the data sample (i.e. fraction of outcomes of the MC 


11 


Abbasi et al.: IceCube Optical Follow-up Program 


experiment which show worse agreement with the background- 
only hypothesis than the actual measurement) is the desired con- 
fidence level for the exclusion of this combination of model pa- 
rameters. 
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